基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。     

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

SSPS-OMEGA超大型球形聚光器结构拓扑构型与设计

为同时兼顾千米量级尺度空间太阳能电站聚光系统的在轨发射安装难度、光学收集效率、稳定性及结构工艺等多方面的因素，提出采用基于球面正多面体的经纬线划分法模块化构建球形聚光器。该方法不仅可实现基础划分单元分布均匀、规格种类少，降低在轨构建与运载发射难度。而且只要合理控制弧高等分数和基础划分单元的口径就可以同时保证聚光系统的光学收集效率。理论与仿真结果显示，当弧高等分数为24时，该划分法共有25种基础单元件，整个球面被分成5762块。并且此时模块化构建的实际球形聚光器的光学收集效率与理想球面相比基本保持不变，解决了太阳能电站聚光系统设计中存在的关键问题。     

整体弹性结构法及在叶片流固耦合分析中的应用

针对三维叶片时域流固耦合振动响应计算普遍耗时的问题,采用一种假设整体结构的模态位移,求解固体时域响应和实现高效网格变形,发展了一种3维时域流固耦合分析的整体弹性结构方法,并应用于压气机叶片0°和180°相位角的气动稳定性分析中。结果表明,所发展方法的计算结果与传统双向时域算法和文献的结果较为吻合,而计算效率相比于传统算法显著提升;在所分析的两个相位角下,叶片振动的气动阻尼均随流量减小先增大后降低,相比于0°相位角,180°下叶片的气动稳定性大幅提高,表明该方法能有效应用于叶片的工程流固耦合研究。     

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

共轴刚性双旋翼/机身干扰流场数值模拟

采用滑移网格技术求解Navier-Stokes (RANS)方程的方法,研究了共轴刚性双旋翼/机身的干扰问题。通过Caradonna-Tung旋翼、Robin直升机、Maryland直升机旋翼/机身干扰和Harrington 2共轴双旋翼等算例,验证了所提出的旋翼流场数值模拟方法的正确性。在此基础上,以Maryland机身为原型,分析了不同桨距的共轴刚性双旋翼与机身之间的干扰特性。结果表明:所提出的数值模拟方法能够很好地模拟共轴刚性双旋翼/机身的气动干扰特性;由于机身对于共轴刚性双旋翼下洗流场的阻滞作用,旋翼的悬停效率增加5%左右,并且随着拉力系数的增大使得悬停效率的增量更加明显;旋翼的悬停效率增加主要来源于下旋翼0°方位角附近的桨叶升力系数的增大,并且拉力系数的增量由桨根向桨尖方向逐渐减小。      

网格拓扑对DLR-F6构型数值模拟的影响

基于TRIP3.0软件平台和多块结构网格技术,开展了控制方程和网格拓扑等因素对DLR-F6构型数值模拟的影响研究。数值计算采用与试验相同的参数,采用了RANS/TLNS方程,生成了O型和H型2种拓扑的粗、中、细网格进行模拟。分别从气动特性、压力系数分布曲线和表面流态3个方面对结果进行分析。通过与试验数据的对比表明,在模拟小分离流动时,采用H型拓扑网格和RANS方程模型,获得的计算结果会更为准确。     

带旋转修正的弹簧-TFI混合动网格方法

网格存在大变形时,传统的超限插值(TFI)动网格方法易于引起网格正交性问题.依据几何关系并考虑其插值特点,提出了一种带旋转修正的TFI动网格方法.对结构网格块进行分块,在子块角点间建立弹簧元,并基于弹簧类比法计算角点位移;子网格内部结点位移采用修正后的TFI动网格方法进行计算,从而形成具有弹簧-TFI混合特征的动网格新...